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We study the stationary probability distribution of a system driven by shot noise. We find that 
both in the overdamped and underdamped regime, the coordinate distribution displays power-law 
singularities in its central part. For sufficiently low rate of noise pulses they correspond to distribu- 
tion peaks. We find the positions of the peaks and the corresponding exponents. In the underdamped 
regime the peak positions are given by a geometric progression. The energy distribution in this case 
also displays multiple peaks with positions given by a geometric progression. Such structure is a 
signature of the shot-noise induced fiuctuations. The analytical results are in excellent agreement 
with numerical simulations. 
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I. INTRODUCTION 

Shot noise is an important source of fluctuations in 
dynamical systems. The discreteness of the modulating 
force which underlies such noise can be due to the quan- 
tization of the electromagnetic radiation that drives the 
system or the discreteness of the electron charge or spin 
in the electric or spin current in the system. The features 
of the system dynamics provide means for identifying the 
noise statistics, which is sensitive to the microscopic na- 
ture of the noise source. An example of using the dynam- 
ics for detecting non-Gaussian noise statistics is the re- 
cent theoretical and experimental work on noise-induced 
switching between coexisting stable states in Josephson 
junctions and mechanical resonators [iHS]. 

The strong effect of the noise discreteness is easy to see 
for a strongly damped noise-driven system localized in a 
potential well. If the noise pulses are all of the same sign, 
the noise pushes the system only in one direction, and the 
stationary probability distribution of the system is equal 
to zero on the one side of its stable state. In the case of 
shot noise, the distribution turns out to be singular near 
the stable state both for the case of pulses with random 
exponentially distributed amplitude flO-13' and pulses 
of constant amplitude 14^. Depending on the pulse rate 
relative to the system relaxation rate, it either displays 
a power-law divergence or goes to zero as a power law of 
the distance to the stable state. 

In this paper we study the probability distribution for a 
shot (Poisson) noise-driven system with arbitrary damp- 
ing. Our results extend from the limit of overdamped dy- 
namics to underdamped dynamics, where the relaxation 
rate is small compared to the typical vibration frequency 
of the system. Examples of underdamped systems that 
are of interest for studying Poisson-noise induced fluctu- 
ations include Josephson junctions 15J, nano-magnetic 
oscillators [1^], and high-Q nanomechanical resonators 
coupled to electron tunneling, to mention but a few; the 
problem attracted much attention recently in the con- 
text of the studies of radiation-pressure shot noise with 
optomechanical systems [T71 [TH]. We assume that the 



Poisson noise pulses have constant amplitude, which is 
relevant for most of the above systems. 

One would expect that, for an underdamped system, 
the aforementioned singularity of the probability distri- 
bution at the stable state should disappear. Indeed, in 
this case the noise makes the system perform random 
vibrations. Therefore, it contrast to an overdamped sys- 
tem, the probability distribution is nonzero on the both 
sides of the stable state. However, we find that the power- 
law singularity at the stable state persists. 

Moreover, we find that the probability distribution of 
Poisson-noise driven underdamped systems can display 
multiple power-law singularities, with positions forming 
a geometric progression. To gain an additional insight 
into this unusual structure we look at the probability 
distribution of the system energy. We show that this 
distribution also has singularities. 

It is somewhat surprising that more is known about the 
tail of the distribution of Poisson-noise driven systems 
[liailllllllllllllillinillO] than about its central part, 
which is singular. Unless the noise is very strong, this 
central part is formed by the motion of the system near 
its stable state, which is generally described by linear 
equations of motion. This allows us to obtain the results 
in a closed form and to study the critical exponents that 
characterize the singularities of the distribution. 

In Sec. [n] we present the model of a simple shot- 
noise driven system with inertia. The onset of multiple 
power-law singularities of the coordinate distribution in 
the overdamped and underdamped regimes is studied in 
Sees. |III| and |IV[ respectively. The onset of the singulari- 
ties of the energy distribution of an underdamped system 
is considered in Sec. [Vj Section VI provides a qualitative 
insight into the onset of singularities and explains their 
positions. Section [VII| contains concluding remarks. 



II. MODEL 

We will consider a standard model where the dynam- 
ics of a system near its stable state is described by the 
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Langevin equation 

q + 2Tq + ijlq = fp{t), 



(1) 



Here, q is the system coordinate counted off from the 
equihbrium position, uJq is the curvature of the effective 
confining potential, and F is the viscous friction coeffi- 
cient. The force fp{t) is a Poisson (shot) noise. It con- 
sists of short pulses which occur at random, the instants 
tn are uncorrelated. The average pulse rate is v. We as- 
sume that fp{t) is independent of the system coordinate, 
in agreement with our intention to keep the leading-order 
terms in q and q. 

The dynamics of the system ([T]) is determined by two 
dimensionless parameters: the relative decay rate F/wq 
and the relative pulse rate v/T. The pulse area g just 
scales the velocity q. For concreteness, we assume g > 0- 

To obtain the stationary distribution of the system we 
assume that the initial state has decayed and write the 
coordinate as 



dt' fp{t')a{t-t'), a{t) 



Ai — A2 



(2) 



In this equation a{t) is the response function of the sys- 
tem. 

Using the well-known expression for the characteristic 
functional of the Poisson noise [H], we can then write 
the stationary distribution of the system coordinate as 

p{q)^{5{q-q{Q)))^ / - exp 



■4){k) = F 



dt 



27r 



ikq- -Tpik) 



(3) 



From Eq. ([s]), ip{—k) = 4'*{k). Therefore, in much of the 
analysis of ^/'(A:) we will focus on the region fc > 0. 



III. POWER-LAW SINGULARITIES OF THE 
COORDINATE DISTRIBUTION: OVERDAMPED 
REGIME 

Integration over the range of small k in Eq. ([3| gives 
a smooth contribution to the distribution p{q). Singular 
behavior of p{q) is determined by the large-fc behavior 
of the function ip{k). We will discuss it in three limiting 
cases, which correspond to the overdamped and under- 
damped regimes, as well as in the critical regime where 
the dynamics changes from over- to underdamped. 

In the overdamped regime the friction coefficient F ex- 
ceeds the frequency ujq. The system does not oscillate in 
the absence of noise. Both eigenvalues Ai,2 in Eq. ^ are 
real. Then one can see that in Eqs. ^ and ([s]) a(t) > 
for t > 0. Given that a{t) — t' for — > 00, function 
^{k) is analytical in the upper halfplane of the complex- 
k plane. For q < one can then add to the integral 



over the fc-axis in Eq. ([3| an integral over a semicirlce 
|fc| 00 in the upper halfplane, which is equal to zero 
for q < 0. Since there are no singularities in the fc-plane, 
the whole integral is equal to zero. Therefore once the 
system becomes overdamped, even in the presence of in- 
ertia p{q) = for g < 0. 

A simple expression for the probability distribution can 
be obtained in the strongly overdamped regime, which 
corresponds to the limit of a large friction coefficient, 
F ^ Wo- In this regime Ai « — 2F and A2 ~ — cjq/SF, 
with I All IA2I; respectively. 



a(t) « (2F)-i [exp (- 



-w2t/2F) 



exp(-2Ft)] 



The main contribution to ipik) comes from the time 
range < ^ 1/F. In this range in Eq. ([3| « 
exp{—LdQt/2T)/2r. One can show that the resulting ex- 
pression for p{q) coincides with the expression for the 
probability distribution obtained in Ref. 14 using a dif- 
ferent method and in the form less convenient for the 
present analysis. We denote function ^p{k) in this ap- 
proximation as V'id'(fc); this function can be expressed in 
terms of the integral cosine and sine. 
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(4) 

where 7e ~ 0.58 is the Euler constant. As seen from 
this equation, parameter k is scaled by the factor T/g. 
Therefore the characteristic width of the distribution (|3| , 
or in other words, the characteristic spatial scale on which 
the system is localized, is ~ ff/F. 
For large fc, from Eq. Q we have 
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(5) 

The large dimensionless parameter used in deriving this 
equation is kg/2T. 

For kg/2V 3> 1 it is important also to keep in ■(/'(^) 

a correction '!/'iV(fc) that comes from the time range 
t < F~^ in Eq. |3|. One can find it by calculating the in- 
tegral over t in E q. ([s] ) for large fc by the steepest descent 



method, see Sec. IV where a similar but more compli- 



cated case is discussed. In the present case the integrand 
has one saddle point, which is located at the extremum 
of a{t) and is given by equation exp(~2Ft) « (wo/2F)^. 
Therefore the corresponding contribution is missed if one 
disregards the term cx exp(— 2Ft) in a{t). The result 
reads 



F /47rF 
Wo V kg 



1/2 



exp 



ikg 

2r 



(6) 



It is seen from this expression and Eq. ([5|) that ijj^^} ex- 
ceeds the fc-dependent correction to the logarithmic term 
in V^od^ for large fc. 
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In the integral over k in Eq. ([3]), we expand 
exp[—(i//T)'iJj{k)] in ip'^J keeping the zeroth and first or- 
der terms. As seen from Eqs. (|3| and ([s]), the zeroth- 
order term in t/'^^'' leads to a power-law behavior of p{q) 
near g = for q > 

p{q)^q-^'"', l3od = l-^ (r»c^o)- (7) 

The exponent /3od is determined by the ratio of the rate 
of Poisson pulses to the relaxation rate of the system 
Wq/F. For Pod < 0, i.e. for sufficiently high pulse rate, 
the distribution p{q) goes to zero for q — > +0. On the 
other hand, for small pulse rate, where /3od > 0, the 
distribution diverges for g — ^ 0. This behavior is seen in 
Fig. [ij which shows the results of numerical simulations 
of the equation of motion (fTl). 



/3od = 1 — (j^/|A2|). The distribution p{q) also has a 
power-law singularity at nonzero q. It is described by 
Eq. ([S]) with the corresponding /3od and with the position 
of the singularity changed from q = g/2T to q = ga{ts), 
where tg is the root of equation Q.{t) = 0. 

The occurrence of the singular peak of the distribution 
for a nonzero g is a benchmark of inertia in overdamped 
shot-noise driven systems. It is clearly seen in Fig. [T] 
For the chosen F/wq — 3 the position of the peak is al- 
ready close to its value .g/2F in the limit of large F/wq. 
We checked that it approaches this value with increasing 
T/luq. The critical exponents of the peaks at q = and 
q ~ 3/2F are in excellent agreement with the analytical 
results. They approach the asymptotic values ([7| and ([s]) 
with increasing T/luq and become within the error of the 
simulations already for T/uiq = 10. 




g/2T Tq/g 

FIG. 1. The results of numerical simulations of the probability 
distribution of a shot-noise driven system in the overdamped 
regime, F/ljo = 3. The data curves 1 to 3 correspond to 
i//r = 0.025, 0.05 and 0.25. The respective values of the crit- 
ical exponent are Pod ~ 0.55,0.1 and -3.5. The inset shows 
the distribution on the logarithmic scale; the slopes of the 
straight lines are given by — /Sod- 

Taking into account in Eq. (jsj) the term oc V'od''' '^^^ 
sees that, for Pod > 1/2, the distribution p{q) has a peak 
at g sa f;/2F. It appears on a smooth background, and 
the difference Sp{q) from the background value is 



We emphasize that the peak is asymmetric: the prefactor 
in 5p{q) depends on the sign of g — .g/2F. The peak ([S]) 
does not emerge if one disregards the inertial term in the 
equation of motion ([T]). For smaller Pod the distribution 
itself does not diverge for q g/^T, but for Pod > — 1/2 
the divergence is displayed by the derivative dgp. 

The analysis can be extended to the case where the 
motion is weakly overdamped: F > wq, but the ratio 
T/uJo is not large. If, as before, < — A2 < — Ai, so 
that a{t) oc exp(A2t) for t — > 00, then to the leading 

order ^'fjik) ~ (F/|A2|)lnfc for k 00. The power- 
law singularity for g — ^ -f is given by Eq. Q, but now 



A. Critical regime 

The results for the overdamped regime can be extended 
and the explicit expressions for the parameters can be 
obtained in the critical regime where the motion changes 
from overdamped to underdamped. This happens where 
|F — wqI ^ r and, consequently, |Ai — A2I ^ |Ai,2| ~ T- 
For F > cjQ the distribution p{q) — for g < 0, whereas 
for T < luq the probability to find the system in the region 
g < is nonzero, but p(g) steeply decays with increasing 
—q for small {ujq — F)/F. 

For F = we have a{t) — texp(— Ft) > 0, and there- 
fore still p{q) = for g < 0. Using the arguments that 
led to Eqs. ([7| and Q one obtains that the power-law 
singularity of p{q) for g +0 has the form q~^" with 
Per = 1 — (I'/T). The distribution p(g) also has a power- 
law singularity of the type of Eq. Q , which is located at 
g = g/eV and is characterized by exponent Per — 1/2. 

For small (wo — F)/F > there emerge additional 
power-law singularities of p(g). However, they are lo- 
cated at exponentially small |g| cx exp[— 7rF/(a;o — F)] 
and therefore are extremely hard to resolve. 



IV. UNDERDAMPED REGIME 

The general expression for the probability distribution 
simplifies also in the case of small relaxation rate, F ^ 
ujQ. In this case the motion of the system in the absence 
of noise is weakly damped vibrations at frequency « wq, 
and A1.2 ~ — F ± itUQ. Noise pulses excite vibrations at 
random. Clearly, the stationary probability distribution 
p{q) in the presence of dissipation and noise is no longer 
limited to the region g > 0, it is expected to be almost 
symmetric with respect to g. 

To describe p{q) we note that in Eq. ^ 

a{t) w ujq^ exp(— Ft) sinwot. 
For an underdamped system the function ip{k) = ipudik) 



in Eq. ([s]) to the leading order in F/c^o has the form 

nkg/ujQ 

^niik) ~ ^nlik) = dx[l-Ux)]/x. (9) 

Jo 

Function -0^,^^ describes the smooth part of the distribu- 
tion p{q). The typical width of the distribution, as seen 
from Eqs. ([S]) and is ~ s/wq. 
In the range kg/tUQ ^ 1 we have 



kg 



(10) 



which indicates that, even in the underdamped case, the 
distribution has a power-law singularity for |g| — > 0. Gen- 
erally, this singularity appears on a smooth background. 
From Eqs. ([S]) and ( 10 ) the difference 5p{q) from the 
background value for small \q\ is 



Spiq)^\q\-^-^, 



/3ud = 1 - 



r 



(r«c^o). (11) 



In contrast to the overdamped case, Eq. Q, the singular 
behavior occurs on the both sides of g = 0, with the same 
exponent /3ud- This exponent is again determined by the 
ratio of the pulse rate u to the relaxation rate, which is 
equal to F for weak damping. 



The approximation leading to Eq. (Ill is justified for 
/3ud > —2. In the opposite case, /3ud < —2, the dis- 
tribution for small q is formed primarily by the region 
where kg/u)Q < 1 and the distribution is parabolic near 
q = 0. Still there may be singularities in the derivatives 
of p of sufficiently high order; in what follows we assume 
/3ud > -2. 



A. Multiple distribution peaks away from the 
origin 

The overall distribution p{q) in the underdamped case 
turns out to have multiple singularities. In their analysis 
one should take into account the terms ^ F/wq in V'ud(fc)- 
In the limit of large k they can be found by calculating 
the integral over time in Eq. ([S]) by the steepest descent 
method. It requires bending the contour of integration 
over time, which is justified since a{t) has no singularities 
near the Re i-axis, see Fig. [2j 

As seen from Fig. [2] for large k the function ipudik) 
has two major contributions. One comes from the re- 
gion of large time, where kga{t) < 1 and the integration 
goes along the Re t-axis. To the leading order in fc, i t 



is given by the logarithmic term in '0ld^(^)i Eq. (10) 



It is not proportional to F/wg ^ 1- The other contrr 
bution comes from smaller times, where the exponential 
term exp[ikga{t)] in il}{k) in Eq. (|3| can significantly dif- 
fer from 1. In this region, if the integration contour is 
appropriately bent, this term has multiple saddle points 
where a{t) = 0. These points are marked in Fig. [2] The 




FIG. 2. The contour of integration over time for calculat- 
ing function V(^)) Eq. ([3|, for an underdamped system for 
kg/ujo 2> 1. In the range exp{~rt')kg /ljo S> 1 (t' = Re t) 
the contour is oscillating and goes through the saddle points 
t„ ^ (n + l/2)7r/tJo of the function exp[ikga{t)], making an- 
gles (-l)"+V/4 with the Re t-axis (n = 0,1,2,...). The 
saddle points are indicated by the solid circles. For large Re t 
the integration contour approaches the Re i-axis. In the over- 
damped limit the integration contour used to obtain Eq. Q 
crosses the Re f-axis only once for exp(— Ft) — ajo/2r and 
then approaches the Re t-axis. 



resulting contribution is oc T/ujq and has the form 

^W(fc)«- 



ujQ \ kg 
</>„(fc) = (-l)"Me-r*._ 

UJQ 

tn = TrujQ\n + l/2). 



n=0 

i-iy 



e""/2exp[z(/.„(A:)] 



(12) 



Here, parameter rtmax is determined by the condition 
[kg/uJo) exp(— Ft„) ^ 1 for n < rimax- Therefore 
IV'id I < 1, and in Eq. ^ ex^[-{v /V)ip\^l] can be ex- 
panded in V'id'; will keep the first-order term in this 
expansion. 

From Eqs. (|3| and ( 12 1 , for sufficiently low rate of noise 
pulses i^/F, where /3ud > 1/2, distribution p{q) has mul- 
tiple power-law peaks. The deviation of p{q) from the 
smooth background near the nth peak 5pn{q) is 



,5p„(g)cx|g-g„r^-+i/2^ 



(-1)' 



(13) 

As seen from Eqs. (12 1 and (131, the positions (?„ of the 
singular peaks of p{q) form a geometric progression. All 
peaks display a power-law shape with the same exponent 
Pud ^1/2. At the same time, the prefactor in 5pn [q) takes 
on different values on the opposite sides of the peak, i.e., 
it depends on the sign of q — Equations ( [l2| and 
( 13 1 make it possible to find the prefactor (this can also 
be done in the overdamped regime), but the expression 
is somewhat cumbersome 

The peaks with largest |g„| are the ones with n = 
and n = 1. They lie on the opposite sides of the center 
of the distribution at g = 0. The positions go and gi 
are asymmetric for nonzero F/wq, but the asymmetry 
is weak for small ttF/ojo- Other peaks lie between go 
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FIG. 3. The results of numerical simulations of the probability 
distribution of a shot-noise driven system in the underdamped 
regime, F/wo = 0.125. The data curves 1 to 3 correspond to 
v/V = 0.15, 0.4, and 0.75. The respective values of the critical 
exponent are /3ud = 0.85, 0.6 and 0.25. The inset shows the 
distribution on the logarithmic scale; the slopes of the straight 
lines are given by /3ud. 
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FIG. 4. The vicinities of the peaks of the distribution p{q) 
that lie away from q = Q for the overdamped (left panel) 
and underdamped (right panel) systems. The distances are 
counted off from the positions of the corresponding peaks g„; 
in the overdamped case there is only one such peak, whereas 
for the underdamped case we chose peaks with n = Q and 
3, Eq. ([13|. The labels "right" and "left" indicate the right 
and left side of the peak, i.e., q — g„ > and g — g„ < 0, 
respectively. In the left panel T /ujq — 10, /3od ~ 0.8; in the 
right panel T/ujo = 0.05, Pud = 0.85. The straight lines show 
the expected asymptotic slopes of In p{q) for q ^ q^. The 
maxima for comparatively large |? — gn| visible in the right 
panel correspond to other peaks of p(g), i.e., the peaks with 
n / 0,3. 



and qi. Their amplitudes increase with decreasing |(7„| 
because of the factor exp(rt„/2). We note again that 
these amplitudes are oc F/wq <C 1. For higher rates of 
shot-noise pulses, where 1/2 > /3 > —1/2, the peaks of 



p{q) disappear in the asymptotic theory (13), but a self- 
similar structure of power-law divergences can still be 



seen in the derivative of the distribution dqp. 

The predicted behavior is in agreement with the results 
of numerical simulations of Eq. (jlj shown in Fig. [S] The 
distribution obtained by simulations displays power-law 
singularities, and the positions of these singularities are 
in excellent agreement with Eq. (13). The exponent of 
the singularity for 5 = is also in agreement with the 
theory. For the peaks of p{q) at g 7^ in the case of 
the smallest v/T, where /3ud = 0.85 and these peaks are 
most pronounced, the exponents agree with the theory. 
However, for the moderately small r/cjg used in Fig. [sj 
for larger v/T the agreement is worse: for /3ud = 0.6 
instead of the expected exponent 0.1, see Eq. (13), the 
obtained exponent is ~ 0.16, which we believe is due to 
the overlapping of the peaks that complicates retrieving 
the exponent. 

Careful studies for smaller F/wq demonstrated an ex- 
cellent agreement of both the positions of the peaks and 
the exponents with the theory. The exponents are the 
same on the both sides of the peaks and are the same for 
all peaks at g„ ^ 0. The results are shown in Fig.|4j This 
figure also demonstrates an excellent agreement with the 
theory of the exponents obtained by numerical simula- 
tions for an overdamped system. 

An important feature of the singularities of the distri- 
bution p{q) is that their positions g„ depend only on the 
area of the noise pulses g and the system parameters, but 
not on the pulse rate u. In contrast, the exponent /3ud 
depends on v scaled by the system relaxation rate, but 
is independent of the noise pulse area g. 



V. SINGULAR PEAKS IN THE ENERGY 
PROBABILITY DISTRIBUTION 

For a Poisson-noise driven system the probability dis- 
tribution over the system energy pe{E) differs from the 
Boltzmann distribution. The tail of this distribution in 
the limit of small damping was discussed earlier |15) . 
Here we also consider the small damping case, F <C wq, 
but we are interested in the central part of the distribu- 
tion. We show that the distribution is singular and can 
have multiple power-law peaks. 

Function pe{E) can be conveniently expressed in terms 
of the distribution of the system in phase space Pq^p{q,p)- 
For a linear system described by the equation of motion 
(nj) the stationary distribution pq^p — {S[q — q{t)]S[p — 
p[t)]) can be found using the explicit expression ([2| for 
the system coordinate q{t) and the corresponding expres- 

sion for p{t) = q{t) = dt' fp{t')dta{t - t'), 
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PEiE) = J dqdpp,^p{q,p)S[E-E{q,p)], E{q,p) = [p" + Lolq^) /2, 



Pq,p{q,p) = (27r) ^ / dkgdkpexp -i^k^q + kpp) - -^^^(k) 



POO 

V'£(k)=r/ dt{l~exp[tgfu{t)]}. 
Jo 



[k — (kq, kp)], 



(14) 



/k(t) = fcga(i) + kpa(t). 



This expression is a straightforward extension of Eq. ^ . 

The singular behavior of the distribution is deter- 
mined by the behavior of the function in the range 



of large jfe^j, \kp\. The leading-order contribution ^/"^^ to 



r 



It is seen from this equation that, for k — >■ cx), the leading- 
order term in ipE is V'^^ ~ 1ii(k(7), which coincides with 
Eq. ( |l5l ). We note, however, that variables (k, 0) are less 
convenient for calculating the behavior of for E 



tpE comes from the region of t in Eq. ( 14 ) where the term than (fcg, kp) because of the singular nature of the in- 

and 



exp[ig/k(t)] is fast oscillating and can be disregarded. 
The size of this region can be estimated by noticing that 
a{t),a{t) (X exp{—rt). Therefore 



JO) 



In 



[.gmax(|fcq|cjo \ |fcp|)] 



(15) 



From Eqs. (14) and (151, the distribution pE has a 
power-law singularity at E ^ 0, 



Pe{E)(xE 



2V 



(E^+O). (16) 



This singularity corresponds to a divergent (but inte- 
grable) peak for Pe > 0; if /3b > —1 the derivative OePe 
diverges for E ^ +0. 



It is interesting to compare the singularity ( 16 ) with 
the singularity of the coordinate distribution p(gj of the 
underdamped system for g — >■ 0, Eq. ( |11[ ). The exponent 
f^E is expressed in terms of the exponent of p{q) as (3e — 
(/3ud + l)/2. Therefore, if /3ud > and p{q) has a peak 
for (7 = 0, the energy distribution also has a peak for 
E — > +0. However, the energy distribution can have a 
peak for E — +0 even where /3ud < and p{q) does not 
have a peak for q = 0. 



A. Multiple peaks of the energy distribution for 
nonzero energy 

The contribution V'^' ipE of the subleading order in 
fcq, kp leads to the onset of singularities in the distribution 
p^; for i? > 0. To find them we first change to cylindrical 
coordinates in k-space. 



^k„ 



KsinO. 



Then, changing in Eq. (14 1 from q,p to the standard 
action-angle variables |22j and integrating over the angle, 
we obtain 



PEiE) = (27r)"^y KdKd9J^{K\/2E)e 



(i^/r)V'E 



(17) 



In variables (k, 6), to first order in F/wq function /k in 
the integrand oi tpE becomes 

/k(t) « Kexp(-Fi) 

X [sin(cLio^ + 0) — (F/wo) sin^sinwoi] . (18) 



tegrals; nevertheless, the above expression for -0^ 
Eq. (17) immediately show that the scaling of p^ for 

E- 



-0 is indeed of the form of Eq. ( 16 ) 



The subleading term ■i/'^^'' can be obtained by calcu- 



lating the integral over time in Eq. ( 14 ) by the steepest 
descent method, as in Sec.|lV) cf. Fig:]2| From Eq. ([l8|, 
to the lowest order in F/wg the saddle points are 



tn{e)^TTUJ^\n+l/2)^UJ-^H 



(19) 



[the instants f„ in Eq. (12) are equal to i„(0)]. 



The result of the integration has the form similar to 
Eq. PI, 



-(FM) {2^lngf^ ^""^'''^^^ni^, 

n 

a„(At, 9) = exp [ighiUe)) - (-l)"i7r/4] . (20) 

Here we have disregarded corrections oc T/ujq unless they 
are multiplied by a large factor. In particular, we keep 
the term oc T/ujq in function /k, since it is multiplied 
by K. It can be large for large k and such tn{9) that 
Kg exp[—Ttn{0)] > loq/T 3> 1. To first order in F/ojo 



MtniO)) 



Ke 



^rt„(o) 



(-1)" 



1 



F 

OJQ 



1 



sin 29 



The condition Kg exp[—Ttn{9)] ^ 1 imposes the upper 
limit on n in the sum over n in Eq. ( |20| . 

To find the most pronounced singularities of pE we 
expand in Eq. (17) exp[— (i^/F)-;/'^^''] to the first order in 
-0^^. Then the calculation of pe reduces to integrating 

t/'^'' over 9, which has to be followed by integration over k 
with the appropriate weight. Since from Eq. ( 14 ) tn{9) > 
0, it is convenient to integrate over 9 from —3tF/2 to tt/2, 
which is seen from Eq. (191 to correspond to n > in 



Eq. (120|). 

For large k, integration of exp[ig fi^(tn{9))] over 9 can 
be done by the stationary phase method. The stationary 
points, dgfk = 0, are located at 9st — mn with integer 
m. At these points 9g/k — 0. Therefore 



d9an{K, 9) 



[Kgyn 



[a„(/t,0)+a;+i(K,0)] , (21) 
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where 



Cn 



{(2r/3c.o)(-l)"e-""(°'} '^'r(l/3)/y3 



(r(x) is the Gamma function). 

Taking into account that Joix) oc x~^/^cos(x — 7r/4) 
for 1, one obtains from Eqs. (17 1 - (21) that pe{E) 



displays power-law singularities for nonzero and near 
an nth singularity the deviation 5pE,n{E) oi pe from the 
smooth background is 

5pEAE)oz\E-E,r^f'-+^/\ 

= ^5'exp[-2n„(0)]. (22) 

As in the case of the singularities of the coordinate distri- 
bution for an underdamped system, for all singularities 
oi Pe the exponents are the same and the positions of 
the singularities form a geometric progression. The ex- 
ponent 2Pe — 4/3 = (2/3) — [v/T) is independent of the 
area of the noise pulses g, whereas the positions of the 
singularities En are independent of the pulse rate v. In 
the most interesting case of comparatively low pulse rate, 
where 2f3E ~ 4/3 > 0, the singularities correspond to the 
peaks of the distribution. We note that the condition 
2f3E — 4/3 > holds for higher pulse rate v than the 
condition /3ud ~ 1/2 > 0, which is necessary for observing 
multiple peaks of the coordinate distribution. 




FIG. 5. The results of numerical simulations of the proba- 
bility distribution over the energy of a shot-noise driven un- 
derdamped system, where F/luq — 0.125. The data curves 1 
to 3 correspond to v/F — 0.15,0.4 and 0.75. The respective 
values of the critical exponent for E — >■ 0-I-, which are shown 
by straight lines in the inset, are /3_b = 0.925,0.8 and 0.625. 

In Fig. [5] we compare the predictions with the results 
of numerical simulations. The positions of the peaks of 
Pe are in excellent agreement with Eq. (22 1 already for 
moderately small T /ujo — 1/8. As expected, the peaks 
become less pronounced with the increasing noise pulse 
rate. The singularities of pE ior E — En > are still 
visible even for 2/3B — 4/3 < (data curve 3 in Fig.|5|; the 
asymptotic theory predicts that, for the corresponding 
/3e, the derivative OePe should diverge for E = En- The 



exponent of the peak for E = +0 obtained numerically 
is also in excellent agreement with the theory. 

In Fig. |6] we present the results of simulations of the 
singularities very close to the peaks of p^;. We find that 
the singularities are well described by the power law, and 
the exponents are in excellent agreement with the ana- 
lytical theory. It should be noted that the peaks found in 
numerical simulations had a structure of doublets, with 
extremely small distance between the peaks in the dou- 
blet. Analytically one might expect a doublet structure 
for not too small T/ujq, as seen from Eq. (21 1. However, 
the corresponding analysis is beyond the accuracy of the 
asymptotic theory developed here. 



^ 10' 
10' 
1 

10-' 



















£3 (right) ^ 



10" 



10" 



10" 



E-E\/s 



FIG. 6. The results of numerical simulations of the probabiHty 
distribution over energy of an underdamped system near the 
peaks ([22| with n = and 3. The labels "right" and "left" 
refer to the right and left sides of the peaks, where E — En > 
and E — En < 0, respectively. The parameters are F/tJo = 
1/20, /3e = 0.925 {u/V = 0.15). The straight line shows the 
analytical value of the exponent of the power-law peaks for 

E^ En. 



VI. QUALITATIVE PICTURE OF THE ONSET 
OF THE DISTRIBUTION SINGULARITIES 

The onset of singularities of the probability distribu- 
tion of Poisson-noise driven systems can be understood 
by noticing that a single noise pulse shifts the momentum 
of the system by 5. If the pulse rate v is small compared 
to the relaxation rate t~^, on average the system has time 
to relax between the pulses and to approach the equilib- 
rium position q = p = 0. This leads to accumulation 
of the probability distribution near the equilibrium posi- 
tion for i>tr <C 1. As we showed, the result is a power-law 
peak of p{q) at q — with exponent 1 — vtr. 

The argument and the result apply to both over- 
damped and underdamped systems; for strongly over- 
damped systems tr = 2r/wQ, whereas for strongly un- 
derdamped systems tr — l/F. The expression 1 — vtr 
for the exponent agrees with the qualitative picture that, 
for — 0, where the pulses are rare, the distribution 
should be proportional to the reciprocal velocity for a 
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given q, which is itself proportional to q (or to \q\, in 
the underdamped case). On the other hand, for large v 
the noise becomes effectively Gaussian and the distribu- 
tion becomes smooth near the maximum, see below. We 
note that even in the large-vtr limit p{q) has a power-law 
cutoff for gq — > +0 in the case of an overdamped system. 

The power-law singularity emerges also in the critical 
regime where the system dynamics changes from over- 
damped to underdamped. A similar singularity, and 
with the same exponent, can be shown to characterize 
the distribution over the momentum p. Interestingly, 
power-law singularities of the distribution emerge also in 
overdamped systems driven by another important type 
of non-Gaussian noise, the telegraph noise [23] . 

An underdamped system excited by a single noise pulse 
performs weakly damped vibrations. If initially the sys- 
tem was at the equilibrium position, the momentum right 
after the pulse is p = 5, and the ensuing vibrations in the 
weak-damping limit have the form 

lit) ~ (3/'^o)e"^* sinwoi (ujq > T); 

here, we count time off from the instant when the pulse 
occurred. 

The slowing down at the turning points q = p = leads 
to the peaks of the probability distribution. From the 
above expression, such peaks should be located at q = 
Qn = (-l)"(g/a;o)exp(-ri„) with t„ = 7rwo('^ + 1/2), 
in agreement with Eqns. (12 1 and (13). As we showed. 



they are described by a power-law with the exponent 
(1/2) — vtj.. This exponent differs from the exponent 
of the peak at q — 0. The onset of the peaks at 7^ 
requires a lower value of vtr than for the q = 0-peak. 

The fact that the exponent approaches 1/2 for 
can be again understood as a result of the probability 
distribution being proportional to the reciprocal velocity 
in this limit. Near the extrema of q{t) the velocity scales 
as \q — q,i\^^^- We note that the peaks of p{q) become 
strongly asymmetric for small vtr, in agreement with the 
above argument. 

In the case of an overdamped system, after the pulse- 
induced increase of the momentum p = — > p = 5, the 
system coordinate first moves away from equilibrium and 
then monotonically comes back. 



q{t) « (.9/2r) (e-o*/2r , 



-2rt 



(^0 «r). 



Respectively, the probability distribution over q has a 
peak where q{t) is maximal, q « g/^T. We showed that 
the peak is described by a power-law with an exponent 
(1/2) - for < 1/2. 

Similar arguments can be applied to the peaks of the 
energy distribution in an underdamped system p^- This 
distribution can have a power-law peak for i5 — >■ due 
to the accumulation of the probability density near q — 
p = for rare noise pulses. To understand the peaks for 
nonzero energies we note that, for low pulse rate, after a 
noise pulse at t = the energy evolves as 

E{t) « Ig^ exp(-2R) [1 - (T/^o) sm2ujot] 



to first order in F/wq- Function E{t) has inflection points 
at tn = nu}Q^{n + 1/2). The peaks of occur as a 
consequence of the slowing down near these points. The 
values En — E{tn) give the positions of the peaks, in 
agreement with Eq. ( 22 ) . The corresponding exponent is 



(2/3) — vtr, the same for all peaks. The limiting value 
of the exponent for z/ — )■ can be again understood by 
noticing that pE in this case is given by the reciprocal 
rate of the change of E[t) near the inflection point. 

In real systems, along with Poisson noise, there are 
present other noises. In particular, if dissipation comes 
from the coupling of the system to a thermal bath, there 
is white Gaussian thermal noise with intensity AVksT, 
where T is the bath temperature. This noise is described 
by an extra force fxit) in the equation of motion ([T]). 
Using the explicit form of the characteristic functional 
for white Gaussian noise, cf. [21], one can show that 
such noise leads to an extra term in the function ^->{k) in 
Eq. (§, 



(i//r)i/;(fc) -> {v/T)ilj{k) + 2kBTujQ^k'^ 



(23) 



In the absence of Poisson noise, v — Q, the distribution 
p{q) ([3| is then Gaussian, with no singularities. The 
distribution is Gaussian also if the noise is Gaussian, but 
not (5-correlated. 



The term oc k"^ in Eq. (23) imposes an effective cut- 



off on the values of k that contribute to the integral 
over k in Eq. (H), |fc| < wo/(fcBr)i/2. Respectively, the 
peaks of the distribution p{q) are described by the power 
laws found in Sees. IIIII and IIVI at distances that exceed 
(ksTy/^ /ujo from the centers of the peaks. This provides 
the bound for the observation of the considered power law 
singularities. 



VII. CONCLUSIONS 

We have studied the central part of the stationary 
probability distribution of shot-noise driven systems. We 
analyzed the probability distributions of the system co- 
ordinate and energy. Our central result is that the distri- 
butions can display singularities, and in particular peaks, 
at the equilibrium position of the system as well as away 
from it. The singularities emerge whether the system 
dynamics is overdamped or underdamped. They are de- 
scribed by the power law, with characteristic exponents. 
The positions of the singularities and the exponents are 
obtained in the explicit form. 

The pattern of the singularities of the stationary dis- 
tribution provides a reliable indication of the presence 
of shot noise. In tunable systems, to observe the singu- 
larities it is advantageous to make the relaxation rate of 
the system tr^ larger than the pulse rate v. Observing 
the peak of the coordinate distribution at the equilibrium 
position requires vtr < Ij whereas the peaks away from 
the equilibrium position emerge for vtr < 1/2. How- 
ever, even where vtr > 1/2, there are still singularities in 
the derivatives of the distribution. The power-law peak 
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in the energy distribution of underdamped systems at 
E = emerges for i/t^ < 2, whereas the peaks for finite 
energies emerge for vtr < 2/3. 

The peaks in the probability distributions away from 
the equilibrium position are characterized by the same 
exponents, which are determined solely by the parameter 
vtr- The exponents are different for the coordinate and 
energy distributions. The positions of the peaks, on the 
other hand, depend on the area of the noise pulses, but 
are independent of the pulse rate v. 

Particularly attractive in terms of identifying the noise 
seem to be underdamped systems, since the distribu- 
tions can display multiple singularities whose locations 
form a geometric progression with common ratio oc 
exp(— ttF/wo). However, if the system is very strongly 
underdamped, the singularities may start overlapping, 
which would complicate resolving them. Also, the char- 
acteristic heights of the distribution peaks decrease with 
the decreasing F/wq, which suggests that T/ojq should 



not be too small. This can be of interest for reveal- 
ing shot noise in side-band cooled vibrational systems, 
since the cooling leads to the increase of the decay rate 

[T71IIH1IM1I2S]. 

The singularities of the probability distribution should 
occur also in more general types of shot-noise driven sys- 
tems than the ones discussed above. An important ex- 
ample of current interest is nonlinear vibrational systems 
modulated by a strong resonant periodic field, with shot 
noise coming from fluctuations of this field or, in the case 
of nanomechanical resonators, charge on the resonator. 
The dynamics of such systems in the rotating frame near 
the stable state of forced vibrations is described by lin- 
earized equations of motion that have a form somewhat 
different from the standard Langevin equation ([T]) [9]. 
However, we expect that the presented results can be ex- 
tended to such systems and to other systems that perform 
small-amplitude fluctuations about their stable states. 
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